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Tensors describing boost-invariant and cylindrically symmetric expansion of a relativistic dissi- 
pative fluid are decomposed in a suitable chosen basis of projection operators. This leads to a 
simple set of scalar equations which determine the fluid behavior. As special examples, we discuss 
the case of the Israel-Stewart theory and the model of highly-anisotropic and strongly-dissipative 
hydrodynamics ADHYDRO. We also introduce the matching conditions between the ADHYDRO 
description suitable for the very early stages of heavy-ion collisions and the Israel-Stewart theory 
applicable for later stages when the system is close to equilibrium. 
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I. INTRODUCTION 

Soft-hadronic observables measured in the ultra- 
relativistic heavy-ion experiments may be very well de- 
scribed by the standard perfect-fluid hydrodynamics (for 
a recent review see [l|) or by dissipative hydrodynam- 
ics with a small viscosity to entropy ratio @-0|- These 
approaches assume generally that the produced system 
reaches a state of local thermal equilibrium within a frac- 
tion of a fermi 1 . 

On the other hand, many microscopic approaches 
assume that the produced system is initially highly 
anisotropic in the momentum space, e.g., see Q. High 
anisotropics present at the early stages of relativistic 
heavy-ion collisions exclude formally the application of 
the perfect-fluid and dissipative hydrodynamics. This 
situation has triggered development of several approaches 
which combine the anisotropic early evolution with a 
later perfect-fluid [9rll2| or viscous [H, [l4| description. 
Very recently, a concise model describing consistently dif- 
ferent stages of heavy-ion collisions has been proposed 
in Refs. (a highly-Anisotropic and strongly- 

Dissipative HYDROdynamics, ADHYDRO), see also a 
similar work that has been presented in Refs. (TqI. [2(jj. 

In this paper, in order to analyze in more detail the 
connections between different effective descriptions of 
very early stages of heavy-ion collisions we consider a 
simplified case of the boost-invariant and cylindrically 
symmetric expansion of matter. In the first part of the 
paper, we introduce tensors that form a suitable basis 
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for decompositions of various tensors characterizing dis- 
sipative fluids. Then, we use this basis to analyze the 
Israel- Stewart [H [22| and ADHYDRO equations [H- 
Il8j ]. Finally, we show how the solutions of the ADHY- 
DRO model may be matched with the solutions of the 
Israel-Stewart theory. The last result may be treated as 
a generalization of the approach presented in where 
no transverse expansion of matter was considered. 

The formal results presented in this paper, when imple- 
mented as numerical procedures, may be used to model 
the behavior of matter produced at the very early stages 
of heavy-ion collisions. Of course, the use of boost- 
invariance and cylindrical symmetry implies that this de- 
scription should be limited at the moment to central col- 
lisions and the central rapidity region. A generalization 
of our framework to more complicated geometries is a 
work in progress. 

II. BOOST-INVARIANT AND 
CYLINDRICALLY SYMMETRIC FLOW 

The space-time coordinates and the four-vector de- 
scribing the hydrodynamic flow arc denoted in the stan- 
dard way as x^ = (t, x, y, z) and 

U» = j(l,v x ,v v ,v z ), 1 = {l-v 2 )- 1 ' 2 . (1) 

For boost-invariant and cylindrically symmetric systems, 
we may use the following parametrization 

U° = cosh 0j_ cosh 77||, 

U 1 = sinh 9± cos <f>, 

U 2 = sinh 9± sin <fi, 

U 3 = cosh 0_l sinh rj\\, (2) 

where 0± is the transverse fluid rapidity defined by the 
formula 

v±_ = \/v 2 + v y = tanh#j_, (3) 
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m is the space-time rapidity, 



'/II 



1, t + z 
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and <f> is the azimuthal angle 



y 

arctan — . 
x 



(4) 



(5) 



In addition to we define three other four- vectors. 
The first one, Z M , defines the longitudinal direction that 
plays a special role due to the initial geometry of the 
collision, 

Z° = sinh?7||, 
Z 1 = 0, 



Z 2 = 0, 
Z 3 = cosh 



(6) 



The second four-vector, X^, defines a transverse direc- 
tion to the beam, 

X° = sinh 9± cosh ^ , 
X 1 = cosh 6*x cos 4>, 



X 



cosh #j_ sin ( 



X 3 = sinh 6*x sinh 7y|| , 



(7) 



while the third four- vector, defines the second trans- 
verse direction 2 , 

y° = o, 

Y 1 = -sin0, 



F = cos< 

y 3 = o. 



(8) 



The four-vector is time-like, while the four-vectors 
Z^jX'-'-jY 11 are space-like. In addition, they are all or- 
thogonal to each other, 



U 2 

u ■ z 
z-x 



1, 



Z 2 = X' 



Y 2 



-1. 



0, U-X = 0, U-Y = 0, 

0, Z ■ Y = 0, x • y = 0. (9) 



All these properties are most easily seen in the local 
rest frame of the fluid element (LRF), where we have 
9±_ = 7711 = <t> = and 

[/ = (1,0,0,0), 

Z = (0,0,0,1), 

X = (0,1,0,0), 

y = (0,0,1,0). (io) 



In the formalism of dissipative hydrodynamics one uses 
the operator A^ v = g^ v — U^U V ', that projects on the 
three-dimensional space orthogonal to U^. It can be 
shown that 3 



= -X^X V - Y^Y V - Z^ l Z v 



(11) 



Using Eqs. ([9]) we find that Z M , X M and Y^ are the eigen- 
vectors of A^, 

A^Xv = x*, A^y = y\ a^z„ = z>*. (12) 



III. EXPANSION AND SHEAR TENSORS 

In this Section we follow the standard definitions of the 
expansion and shear tensors [13] and show that they can 
be conveniently decomposed in the basis of the tensors 
obtained as products of the four-vectors X 11 , Y^ and Z*. 

The expansion tensor is defined by the expression 



V = A£A^£/ Q) , 



(13) 



where the brackets denote the symmetric part of dpU a . 
Using Eqs. ([2]) in the definition of the expansion tensor 
(fl~3|) and also using Eqs. ©-((U, we may verify that the 
following decomposition holds 



6» v =6 X X»X» +6 Y Y»Y V +6 z Z* i Z v , (14) 



where 



ci9 f!9 
9 X = X^X^" = — ^ cosh(9 ± - sinh(9 ± , (15) 



Or 



8y = YaY v 6^ 



Ot 



sinh 6 i 



and 



7 Z = Z»Z V 8» V = - 



cosh 9 i 



(16) 



(17) 



The contraction of the tensors A^ and 9^ v gives the 
volume expansion parameter 



9 = A^9^. 
Equations CET]) and ([T5J yield 



1X — VY — VZ- 



(18) 



(19) 



Substituting Eqs. fI5 ]l -([T7 | in Eq. (JTHJ) we find that this 
formula is consistent with the definition 9 = d^U^. 
The shear tensor is defined by the formula 

(20) 



7 [iij @tiv ^ApuQ. 



2 By analogy to the terminology used in interferometry (HBT) 
studies, one may say that the four-vector defines the out 3 One may check easily that holds in LRF. Hence, as a tensor 
direction, while Y M defines the side direction. equation, M f II should hold in all reference frames. 
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With the help of the decompositions (fTTj) and (IT4l) we 
may write 



= axX^X" + ayY^Y" + ozZ^Z", (21) 



where 



3 3r 



cosh 9± sinh 9± 



3r 



(22) 



289 ± . 2d0 ± 
--— — smh6»i - -— — cosh6> ± , 
3 or 3 or 



cosh 9 1 2 sinh ( 



3 3r 



3r 



(23) 



and 



+ -— sinhfli + -— — cosh 6* j 
3 or 3 or 



2 cosh 9± sinh#^ . „, 
>z = s — - + (24) 



3r 



3r 



+ ^ sinh 6± + - — — cosh ± . 
3 or 3 Or 

In agreement with general requirements we find that 

ox + o Y + o z = 0. (25) 

In the case where the radial flow is absent 
ox = °Y = 1/(3t) and oz = — 2/(3r), which agrees 
with earlier findings [22j. 

IV. ENERGY-MOMENTUM TENSOR 

The energy-momentum tensors of the systems consid- 
ered below in this paper have the following structure 

= eU^U" + P X X^X U + PyY^Y v + P Z Z^Z V '.(26) 

The quantity e is the energy density, while Px , Py and 
Pz are three different pressure components. In LRF the 
energy-momentum tensor has the diagonal structure, 



e 



Py 

P z 



(27) 



Since we consider boost-invariant and cylindrically sym- 
metric systems, e,Px,Py and Pz may depend only on 
the (longitudinal) proper time 



Using the form of the energy-momentum tensor (|26|) in 
([3"Uf and projecting the result on U u , Z u , X v and Y v one 
gets four equations 

e + ed^U^ + P z U v Z^d^Z u (31) 

+P X U V X^X V + PyllvY^Y" = 0, 
eZ v U v - Z^d^Pz - PzdpZi* (32) 

+P x Z v X»d, i X v + PyZ v Y^d^Y v = 0, 
eXj] v - X^d^Px - Pxd^ (33) 

+P Y X U Y^Y V + p z x v z»d^z v = 0, 
eY v U v - Y^d^Py - Pyd^Y^ (34) 

+P x Y v X fi d tl X v + P z Y u Z^d^Z v = 0. 

Here the dot denotes the total time derivative (the oper- 
ator U a d a ). With the help of the relations 



Y»d^ = 



d^ L = 0, z v = 0, 



dfjY^ = 0, Y v = 0, 



X"d u Z u =0, X»d, l Y v =0, 



(35) 



one can show that Eqs. (|3"!?)) and (IMl) are automatically 
fulfilled. In this way, we are left with only two indepen- 
dent equations: 



^cosh 6»_l^- 
+e cosh 9± 
+P X I coshfl 



sinh 9± ^- I : 

Or , 



(36) 



1 89±_ 
t dr 
d9 A 
dr 



sinh 9 



sinh 9j_ 
88 ± 



89 ± 



sinh 9 1 cosh 9 1 
Py -+Pz 







and 



a 9 

sinh 0j_ — + cosh 9 ± — ) Px 
or or , 



(37) 



(89 
sinh 8j_ — — 
8r 



cosh 9 



-Px 
-Py 



sinh 9± 
cosh 9 



1 00j_ 

t cV 

sinh 
Pz 



dO_L 

4- cosh #j 
= 0. 



«96> ± 



(28) 



V. ISRAEL-STEWART THEORY 



and radial distance 



= \J x 1 + y 2 . 



(29) 



The hydrodynamic equations include the energy- 
momentum conservation law 



(30) 



A. Stress tensor 

In the Israel-Stewart theory , the crucial role is played 
by the stress tensor tt^ that satisfies the following dif- 
ferential equation [U [22| 



(38) 
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Here t„ is the relaxation time, r\ is the shear viscosity, 
and F n is our abbreviation for the scalar quantity 



(39) 



where T is the temperature and a,\ is one of the kinetic 
coefficients appearing in the Israel-Stewart theory. 

The structure of the shear tensor, Eq. (PJTj) . suggests 
that we may seek the solutions of Eq. (155)) in the form 
analogous to Eqs. (|T4l) and (|21j) . namely 



= -kxX^X" + ■k y YI 1 Y v + Tt Z Z fl Z U . (40) 



The condition A^tt^ — leads to the constraint 



TTx + 7Ty + TTZ = 0. 



(41) 



The time derivative of tt^ generates nine terms. Since 
X^ L ,Y^ and Z M are the eigenvectors of the projection 
operator A^ v , see Eq. (fl"2)) . we find 

A£A^ a/9 = 

tt'xX^Xu + 7r x A^X Q X„ + uxX^A^Xp 

(42) 

For boost-invariant and cylindrically symmetric systems, 
the explicit calculations show that the terms A^X a , 
A"y a , and A^Z a vanish. Therefore, the decomposition 
(|4*0)l is indeed appropriate, and Eq. (j3"8")l splits into three 
scalar equations 



Trr^X + TTX = + F V TTX, 

T-nTTy + TT Y = 2f]<7 Y + F^TTy , 

TttTTz + ~nz = Zrjaz + F^ttz- 



(43) 



Due to the constraints ([25]) and (|44j) . only two equations 
in (j4"3")l are independent. 



B. Bulk viscosity 

The isotropic correction to pressure, IT, satisfies the 
equation 



r n n + IT = -(6 + F C U, 



(44) 



where rn is the relaxation time for II, C is the bulk vis- 
cosity, and 



(45) 



where is another kinetic coefficient appearing in the 
Israel-Stewart theory. 



C. Energy- momentum tensor and hydrodynamic 
equations 

In the Israel-Stewart theory, the energy momentum 
tensor has the following form 



(46) 



where P cq is the equilibrium pressure connected with the 
energy density by the equation of state, P eq = P eq (e), 
and II is the isotropic correction to pressure. A simple 
comparison of Eqs. (|46|) and (|26|) leads to the identifica- 
tions 

P X = P cq + U + TTx,P Y =Poq+n + 7Ty, 

P Z = Peq + n + 7T Z = P eq + II - TT X ~ TTy. (47) 



Substituting Eqs. gTJ into Eqs. (H and §7$ we 
obtain two equations for five unknown functions: e, 9x, 
II, ttx, and tty- The two first equations in (|43p as well as 
Eq. (I44p should be included as the three extra equations 
needed to close this system. 



VI. ADHYDRO MODEL 

The Israel-Stewart theory describes the system that is 
close to local equilibrium. Formally, this means that the 
corrections to pressure (the quantities II, ttx, and tty) 
should be small compared to P q- This condition cannot 
be fulfilled at the very early stages of heavy ion collisions. 
In the limit r — > the components of the shear tensor, see 
Eqs. (f2"2"l) - ([2"4"]) . diverge and induce very large changes of 
TOCi tty, and ttz through Eqs. (14"3"1) . This leads to strong 
deviations from local equilibrium. 

In this situation, one tries to construct phenomenolog- 
ical models of the very early stages which grasp the es- 
sential features of the produced matter and may describe 
effectively the early dynamics. On the basis of micro- 
scopic models of heavy-ion collisions, we expect that the 
system formed at the very early stages of heavy-ion colli- 
sions is highly anisotropic — the two transverse pressures 
are equal and much larger than the longitudinal pressure. 

Such anisotropic systems are described most often by 
the anisotropic distribution functions which have the 
form of the squeezed or stretched Boltzmann equilibrium 
distributions for massless partons (this is often called the 
Romatschke-Strickland ansatz (23|) 



gexp 



-JLy/(p.U)* + (x-l)(p-Z)< 
A i 



.(48) 



Here p is the particle's four momentum, g denotes the 
number of the internal degrees of freedom, Aj_ may be in- 
terpreted as the temperature of the transverse degrees of 
freedom, and x is the anisotropy parameter. The energy- 
momentum tensor of the system described by the distri- 
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bution function (gHJ has the form [IMLm 4 
T" u = (e + P±) U"U U - P± gT - (P ± - P^Z^Z". (49) 
This form agrees with Eq. (|26|) if we set 

Px=Py = P±, Pz=P\\- (50) 



In the case described by the distribution function 
the energy density and the two pressures may be ex- 
pressed as functions of the non-equilibrium entropy den- 
sity a and the anisotropy parameter x [Tc| , 



e(a,x) = e eq (a)r(x), 



(51) 



Px(tr,x) = P e(l {a)[r{x) + 3xr'{x)}, (52) 



P||(a,z) = P eq (cr) [r(x) - 6xr'{x)] 



(53) 



We emphasize that s C q(<j) and P C q(cr) are equilibrium 
expressions for the energy density and pressure [l| but 
the argument is a non-equilibrium value of the entropy 
density 



t \ nr. / \ 3fl / 7T 2 (T 
£ eq (a) = 3P eq (cr) = ^ ( 



The function has the form 

x arctan — 1 



r(x) 



x 3 



1 



(54) 



(55) 



Substituting Eqs. <|5D j> — (j55 ]> into Eqs. <(36j) and (|57|l 
(which follow directly from the energy-momentum con- 
servation law) we obtain two equations for three unknown 
functions: cr, x, and 9±. The third equation follows from 
the ansatz describing the entropy production in the sys- 
tem, 



The entropy source £ is taken in the form 



£(<7, x) 



(56) 



(57) 



' cq 



The quantity r cq is a timescale parameter. The form (|57|) 
guarantees that £ > and £(<r, x = 1) = 0. We stress 
that ADHYDRO accepts other reasonable definitions of 
the entropy source. In particular, it would be interest- 
ing in this context to use the forms motivated by the 
AdS/CFT correspondence [Hj]. 



In the original papers the four- vector is denoted as W*. 



VII. MATCHING CONDITIONS BETWEEN 
ADHYDRO MODEL AND ISRAEL-STEWART 
THEORY 

In this Section, we show how the initial evolution of 
the system described by the ADHYDRO model may be 
matched to a later non-equilibrium evolution governed 
by the Israel-Stewart theory. 

We propose to do the matching at the transition proper 
time Tt r when the anisotropy parameter x becomes close 
to unity in the whole space, i.e., when the condition 
\x(t = rt r , r) — 1| -C 1 is satisfied for all values of r. Our 
earlier calculations done within the ADHYDRO frame- 
work [l8| show that if the initial value of x is independent 
of r, the later values of x depend weakly on r, hence, it 
makes sense to use the value of the transition time that 
is to large extent independent of r. Certainly, one should 
always check the sensitivity of the obtained results with 
respect to the chosen value of r tr . The acceptable results 
should exhibit weak dependence on T tr . 

We emphasize that the proposed matching procedure 
differs from our previous strategy where the ADHY- 
DRO model was used to describe the whole evolution 
of the system; from a highly- anisotropic initial stage 
to hadronic freeze-out p^H18j. The use of the ADHY- 
DRO model alone implies a smooth switching from a 
highly- anisotropic phase (where x ^> 1 or x <C 1) to 
the phase described by the perfect-fluid hydrodynamics 
(where x m 1). Moreover, as it has been shown in [l5l[l9| . 
for purely longitudinal and boost-invariant expansion 
of matter, ADHYDRO agrees with the Israel- Stewart 
framework in the intermediate region where \x — 1| <C 1. 
On the other hand, if the transverse expansion is in- 
cluded, the presence of non-negligible shear viscosity trig- 
gers differences between the two components of the trans- 
verse pressure, Px and Py . This effect is not included 
in the ADHYDRO model. Therefore, in the approaches 
that include transverse expansion and noticeable effects 
of viscosity, the matching proposed below is in our opin- 
ion more appropriate than the use of ADHYDRO alone. 
However, the use of the ADHYDRO model is reasonable 
for the situations where the effects of viscosity at the later 
stages of the collisions may be neglected. 



A. Energy-momentum matching 

In order to connect the solutions of the ADHYDRO 
model with the solutions of the Israel-Stewart theory, we 
demand that all thermodynamics- and hydrodynamics- 
like quantities are continuous across the transition 
boundary fixed at the transition time Tt r , namely: 

i) the energy density e is the same on both sides of 
the transition and identified with the equilibrium energy 
density in the Israel-Stewart theory, e — e cq , 

ii) the transverse flow, quantified by the value of 8± , is 
the same at the end of the ADHYDRO stage and at the 
beginning of the Israel- Stewart stage, 
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iii) also the three components of pressure are the same, 
namely 



P 



n + Tix = p± 

Peq + n + TTy = P ± 



eq 



(58) 



P, 



n + 7r z = p. 



Here the values on the right-hand-side are obtained at 
the end of the ADHYDRO evolution, at r = r tr , and 
treated as the input for the stage described by the Israel- 
Stewart equations for r > r tr . From (1581) we find first 
that 7Ty = 7r x and nz = —2t:x- In the next step we find 
that 



KX 



£ C q(CT)xr'(a;) W £ oq (cr eq ) 



4(a; - 1) 



45 



and 



Peq + n 



2Pi +Pn 



3 



(59) 



(60) 



For massless particles considered here, P cq = e cq /3 which 
implies II = 0. 

Finally, we demand that iv) the entropy density is the 
same before and after the transition. This is discussed in 
more detail in the next Section. 



B. Entropy matching 

It is interesting to see in more detail how the last con- 
dition is realized in practice. First, we start with the 
ADHYDRO formulation. If the energy density e corre- 
sponds to the equilibrium energy density, see the con- 
dition i), the corresponding equilibrium entropy density 
may be obtained with the help of the inverse function to 
the function e cq (cr). In this way, expanding the function 
r(x) at x — 1, we find 



J oq — £ cq 



e eq (a)r(x)} 
. q (a) + e eq (cr) 



2(x — I) 2 



c + [de cq (a)/da] s eq (a) 



45 

2(x- l) 2 



(61) 



45 



In the leading order in deviations from the equilibrium, 
we may replace a by a cq in the second term in the last 
line of (1611) . Using the thermodynamic identity for the 
system of massless particles we find 



[de cq (a cq )/da cq ] 1 e cq (o- cq ) 



l^eq) 



cqv^cqy/ u,L/ cqJ °cq\ iy cq/ ^ 

Combing the last two results we find 

{x-lf 



3 

^eq- 



* eq 



l 



30 



(62) 



(63) 



Now we calculate the connection between the equi- 
librium and non-equilibrium entropy density using the 
Israel- Stewart theory. The basic relation in this context 
has the form Wa 



n 2 -k Xv -k 
a = (T cq - fi — - P2— 



A// 



(64) 



In our case we have /?o = Ti/C: fi% = r K /(2rj), and 
7TAi/T = 67r^. We also use the result 



5?/ 
Ta eq ' 



(65) 



which has been derived in Ref. [13|. Following [13j we 
note that expansion of anisotropic distributions around 
the equilibrium backgrounds does not lead to the situa- 
tion described by the 14 Grad's ansatz, hence (1551 differs 
from the standard result by a factor of 6/5. 
The relations listed above allow us to write 



a = a 



eq 



1 - 



135 tt| 
32e 2 q 



r n H 2 



2CT 



(66) 



Using Eq. (|59p in this equation we find 



1 - 



30 



r n H 2 
2CT 



(67) 



Thus, we conclude that matching between ADHYDRO 
and the Israel-Stewart theory is continuous if the condi- 
tion (|65|) is fulfilled and n = at the transition time (in 
agreement with our remarks following Eq. (1601) ). 

We note that the condition (|65[) guarantees that the 
entropy production in ADHYDRO has the same form 
as in the Israel-Stewart theory, which has been already 
shown in Refs. [l5|,ll9|. We also note that for r > r tr , the 
two transverse pressure start to differ from each other, 
since their dynamics is governed by different components 
of the shear tensor. 



VIII. CONCLUSIONS 

In this paper we have introduced a basis of projection 
operators which allows for simple analysis of dissipative 
fluid dynamics of boost-invariant and cylindrically sym- 
metric systems. We have used this basis to analyze the 
equations of the Israel-Stewart theory and the ADHY- 
DRO model. We have shown how the very early evolu- 
tion of matter produced in heavy-ion collisions may be 
described by the ADHYDRO equations combined with a 
later Israel-Stewart dynamics. 
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